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Table 1: List of ensembles used in this study, with uq the tadpole factor and r\ja the scale from the heavy 
quark potential. The r\ja values shown come from a smooth interpolation. 

1. Introduction 

The MILC collaboration has been carrying out simulations of 2+1 flavor lattice QCD with 
an improved staggered quark action for about 10 years. The physics program has recently been 
reviewed in Ref. [jl|]. An important aspect of the MILC collaboration's research program has been 
the study of the light pseudoscalar meson sector. Here we give the latest update of this program. 
Compared to the last status report in Ref. [||] lattice ensembles with smaller lattice spacings, smaller 
light quark masses and lighter-than-physical strange quark masses are analyzed. Furthermore, we 
do fits based on both SU(2) and SU(3) chiral perturbation theory (ZPT), rather than just SU(3) 
as before, and we now include NNLO chiral logarithms. The SU(2) chiral fits are described in 
Ref. %. 

2. The ensembles and the fitting procedures 

The MILC collaboration has generated lattice configuration ensembles at six different lattice 
spacings, ranging from a sa 0. 1 8 fm down to a fx 0.045 fm. In the present analysis, only the a pa 0.09 
fm ("fine"), a w 0.06 fm ("superfine") and a sa 0.045 fm ("ultrafine") ensembles are considered. 
With our very precise numerical data, adding in coarser lattice spacings would require inclusion of 
higher order discretization effects in the fits, which is currently not feasible. 

The ensembles considered in this study are listed in Table [j]. In our notation, am' is the sim- 
ulation light quark mass, with up and down quark masses being equal, and am' s is the simulation 
strange quark mass. Notice that several ensembles have an unphysically light am' s , about 60% 
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Table 2: Masses (in MeV, using r\ = 0.31 17 fm for the scale) for the lightest sea-quark pions of various 
tastes at each lattice spacing. The Goldstone pion is the taste pseudoscalar and has the lightest mass of all 
tastes, while the taste singlet has the heaviest mass. The root -mean-squared (RMS) mass is the average that 
is used in the NNLO chiral logarithms. Unless otherwise indicated, the masses given are also the lightest 
valence-quark pions on each ensemble at that lattice spacing. We drop the a w 0.15 fm and a » 0.12 fm 
ensembles from the current analysis because of the large splittings and heavy singlet pions. 

of the physical strange quark mass, and one ensemble has three degenerate (light) quarks. These 
ensembles were created specifically to have good control over the SU(3) ZPT fits. 

We determine the scale r\ on every ensemble from the static quark potential (see Ref. [[I]]). 
The values listed in Table [| come from a smooth interpolation. For the analysis presented here, 
however, we use a mass independent scheme, where r\ is taken from the smooth interpolation with 
the quark masses set to their physical values. This procedure avoids spurious dependence on the 
quark masses in the ZPT fits. 

Even with the use of the improved staggered (asqtad) fermions and the fairly small lattice 
spacings considered, the taste-violation lattice artifacts are significant, and need to be accounted 
for in the analysis. We do this, as in our previous studies, by using rooted staggered ZPT forms 
(rSZPT) at NLO in our chiral fits [Q, ||]. The "rooting procedure," taking the fourth root of the 
fermion determinant when generating the lattices, is used to eliminate the unwanted tastes present 
with the use of staggered fermions. As reviewed in Ref. [jl|], recent work suggests strongly that the 
procedure does indeed produce the desired theory in the continuum limit. 

As a new feature in the present analysis, our XVI fits now include the NNLO chiral logarithms 
derived by Bijnens, Danielsson and Lahde [Q, ||]. In contrast to the NLO chiral logs, however, 
lattice artifacts are not included in the NNLO chiral logs. Instead, we use the root mean square 
average (over tastes) pion mass for the argument of the NNLO chiral logs. This is systematic at 
this order in ZPT only if chiral symmetry violations from taste-violating lattice effects are signif- 
icantly smaller than the usual chiral violations from mass terms. That begins to be true for the 
a as 0.09 fm points, and is better satisfied for the a 0.06 and 0.045 fm ensembles. It is not true 
for ensembles with a > 0.12 fm, which is why that data is omitted from the analysis. Table ||] gives 
some representative pion masses for our ensembles. 

The SU(3) chiral fits are done in two stages. The first consists of "low-mass" fits used to 
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Figure 1: Low-mass SU(3) chiral fits. The red line is the continuum limit with (light) valence and sea quark 
masses set equal and the strange quark mass fixed at 0.6m% ys . 




0.00 0.05 0.10 0.00 0.05 0.10 

(m.+m^r, x (Z m /Z» ne ) (m.+m^r, x (Z m /Z» ne ) 

Figure 2: Test of convergence of SU(3) XYY fits in the continuum, with the strange quark mass fixed at 
0.6mf ys . For this test we also included NNNLO analytic terms in the fit. 

determine the LO and NLO low energy constants (LECs), namely what we call fa and Bt, (at LO) 
and the Gasser-Leutwyler parameters L, (at NLO). Here the goal is to keep only those ensembles 
and valence points where meson masses (including kaons, which have a quark of mass m' s ) are 
sufficiently light that SU(3) ZPT may be expected to be rapidly convergent. In addition, taste 
splitting as a fraction of the Goldstone pion mass should be small enough that omission of taste- 
violations from the NNLO terms (but inclusion at NLO) is systematic; this, for example, is another 
reason to drop the a ps 0.09 fm ensemble with am' = 0.00155, am' s = 0.031. After these cuts, only 
the three fine and one superfine ensembles with m' s < 0.6mf /!) S are included, and the valence masses 
are limited by m x + m y < 0.6m$ hys . The fits are illustrated in Fig. [j] To test convergence, the full set 
of NNNLO analytic terms may also be added; as shown in Fig. |2[ the convergence is satisfactory. 
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Addition of such terms does not improve the goodness of fit, as can be seen by comparing the 
confidence levels (CL) of the two fits in Figs. [I] and ||[ The fits include all partially quenched data 
for pion and "kaon" (with lighter than physical strange quark mass) decay constants and masses. 

In the second stage, the "high-mass" SU(3) XVT fits, all ensembles listed in Table [I] are in- 
cluded with the valence masses restricted to m x + m y < \.2m phys . The LO and NLO LECs are fixed 
at the values from the low-mass fits. NNNLO and NNNNLO analytic terms are included, but not 
the corresponding logs. These terms are needed to obtain good confidence levels, and they allow 
us to interpolate around the (physical) strange quark mass. The fact that they are required indicates 
that SU(3) ZPT is not converging rapidly at these mass values, unlike the situation in the low-mass 
case. Since the LO and NLO LECs dominate the chiral extrapolation to the physical point, the 
results for decay constants and masses are insensitive to the form of these NNNLO and NNNNLO 
interpolating terms, as long as the fits are good. The high-mass fits are used to give the central 
values of the physical decay constants and other quantities involving the strange quark mass, such 
as fz, B 2 and chiral condensate (uu)2, which are defined in the two-flavor chiral limit (m — ► 0, m s 
fixed at m p s hys ). The high-mass fits are illustrated in Fig. |[ 

3. Preliminary results 

In a first analysis we use, as before, a lattice scale determined from T-splittings ^ which leads 
to rf ys = 0.318(7) fm With this, we obtain 



f„ = 128.0 ± 0.3 ± 2.9 MeV , 

f K = 153.8 ±0.3 ±3.9 MeV , (3.1) 
fx/ fx = L201 (2) (9) . 

Here, and in the following results, the first error is statistical and the second is systematic. 

Our result for f n agrees nicely with the latest PDG 2008 value, f % = 130.4 ± 0.2 MeV 
Since f K is our most accurately determined dimensionful quantity, we can use it to determine the 
scale. This gives r^ hys = 0.3 117 (6) (j^j) fm. Redoing our analysis with this more accurate scale, 
we obtain 

f K = 156.2 ± 0.3 ± 1.1 MeV , f K /f n = 1. 198(2) (If ) , 

f 2 = 122.8 ±0.3 ±0.5 MeV , B 2 = 2.87(1) (4) (14) GeV , 

f 3 = 110.8 ±2.0 ±4.1 MeV , fl 3 = 2.39(8) (10) (12) GeV , 

Ulh = 1-062(1)(3) , hlh = 1.172(3)(43) , 

(uu) 2 = -(279(1) (2) (4) MeV) 3 , (uu) 3 = -(245 (5) (4) (4) MeV) 3 , 

2L 6 -L 4 = 0.16(12)(2) , 2L 8 -L 5 = -0.48(8) (21) , (3.2) 

U = 0.31(13)(4) , L 5 = 1.65(12) (36) , 

L 6 = 0.23(10)(3) , L 8 = 0.58(5) (7) , 

m s = 89.0(0.2)(1.6)(4.5)(0.1) MeV , m = 3.25(1) (7) (16) (0) MeV , 

m u = 1.96(0)(6)(10)(12) MeV , m d = 4.53(1) (8) (23) (12) MeV , 

m s /m = 27.41 (5) (22) (0) (4) , m u /m d = 0.432(1) (9) (0) (39) . 
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Figure 3: High-mass SU(3) chiral fits: in the top plots selected partially quenched data points are shown, 
while in the bottom plot only full QCD points, i.e., points with sea and valence quark masses set equal, are 
shown. 



Here the NLO LECs L, are in units of 10~ 3 , evaluated at chiral scale m r] , and the LO LECs Bj, 
quark masses and chiral condensates are in the MS scheme at 2 GeV. For the conversion from the 
bare quantities we use the two-loop renormalization factor of Ref. [12]. The resulting perturbative 
error is listed as the third error in these quantities. The subscripts "2" and "3" refer to the two-flavor 
(with m s at its physical value) and three-flavor chiral limits, respectively. The quark condensates are 
related to the LO LECs by (uu) j = —fjBj/2. Quark masses, finally, have a fourth error, accounting 
for our limited knowledge of electromagnetic effects on pion and kaon masses (see Ref. [13] for 
how we address this). 

We note that our new results for the decay constants, quark masses, and condensates agree, 
well within errors, with our previous analysis using NLO SU(3) XVI supplemented by higher-order 
analytic terms [0]. Most also have smaller errors. Not surprisingly, however, some of the NLO 



6 



Results from MILC's SU(3) chiral perturbation theory analysis 



Urs M. Heller 



LECs changed considerably with the inclusion of NNLO chiral logs. Similar changes have been 



observed in continuum extractions of these NLO LECs; see for example Ref. [|14fl. The comparison 
with our previous results suggests that NLO SU(3) ZPT plus analytic terms, when implemented in 
a careful manner, can be used to reliably extrapolate physical quantities such as light pseudoscalar 
meson decay constants, , and heavy-light meson decay constants and form factors to the physical 
light quark masses and continuum [jj]]. 

From the ratio of fx/ fn in Eq. we can obtain 

|\g =0-2247(11!) , (3.3) 

which is a significant improvement over our previous result, \V US \ = 0.2246(lj3) 



Using one-loop conversion formulae J15Q we obtain from the SU(3) NLO LECs in Eq. ( p.2| ) 
the scale invariant SU(2) NLO LECs ^ 

h =3.32(64) (45) , k = 4.03(16)(17) . (3.4) 

We observe nice agreement between the SU(3) chiral fit results described here and the results 
of the SU(2) chiral fits given in Ref. [J3j] for all quantities that can be directly compared, namely f n , 
f 2 , B 2 , m, (uu) 2 and J 3)4 . 
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